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Amplitude death can occur in chaotic dynamical systems with time-delay coupling, similar to the 
case of coupled limit cycles. The coupling leads to stabilization of fixed points of the subsystems. 
This phenomenon is quite general, and occurs for identical as well as nonidentical coupled chaotic 
systems. Using the Lorenz and Rossler chaotic oscillators to construct representative systems, 
various possible transitions from chaotic dynamics to fixed points are discussed. 
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I. INTRODUCTION 

Amplitude death occurs when coupled oscillators drive 
each other to fixed points and stop oscillating 0, 0, Q • 
Initial studies of this phenomenon were on systems 
wherein the interaction between the subsystems were as- 
sumed instantaneous. In the absence of coupling the dy- 
namics in the individual subsystems were on limit cycles 
of different periods, near Hopf bifurcations. New stable 
fixed points are created by the coupling, and these be- 
come the attractors of the dynamics. 

There has been considerable recent interest |2[ on the 
effects of coupling in a variety of nonlinear systems, 
not only in the physical sciences, but also in biologi- 
cal and social systems in which in addition to the possi- 
bility of chaotic dynamics, complex phenomena such as 
synchronization, spatiotemporal intermittency, hystere- 
sis etc. have been investigated. Theoretical as well as 
experimental studies of coupled systems have addressed 
a number of issues which include phase-shifting, phase 
locking as well as amplitude death. 

It is often necessary to take account of inherent time 
delays in the coupling, and Reddy et al. [j| investigated 
the collective dynamical behavior of limit-cycle oscilla- 
tors interacting diffusively through time delay at a Hopf 
bifurcation (see also Sec. Illlf) . They observed amplitude 
death of oscillations, regardless of the frequencies of indi- 
vidual oscillators, both theoretically and experimentally. 
In other work, amplitude death phenomenon has been 
studied theoretically in distributed time-delay [(| and in 
a ring of delayed-coupled Q oscillators and experimen- 
tally in coupled electronic circuit [8| and coupled thermo- 
optical oscillators 0- Phenomena similar to amplitude 
death are also observed in time-delayed self-feedback in 
limit-cycle oscillators ^(j and chaotic oscillators [Tl| . 

That this phenomenon can occur for coupled chaotic 
oscillators as well is the main result of the present work. 
By introducing time-delay in the coupling of chaotic sub- 
systems, irrespective of their type, one can achieve am- 
plitude death. This transition can be direct, namely from 
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chaotic motion to fixed points, or indirect, when first a 
pair of limit cycles (rather than a pair of fixed points) are 
stabilized through the coupling. The motion then asymp- 
totically becomes periodic, and the subsystems oscillate 
at a common frequency. As a function of the time-delay, 
however, this limit cycle motion can go from being in- 
phase to being out-of-phase. This transition is signaled 
by a dramatic change in two indicators: the phase dif- 
ference between oscillators, and their common frequency. 
Further variation of the time-delay eventually leads to 
amplitude death via stabilization of new fixed points. 

In the following section we show that identical coupled 
chaotic systems with matching or mismatched param- 
eters, and indeed even completely different chaotic sys- 
tems, all show amplitude death. The in-phase to out-of- 
phase transition in coupled limit cycle oscillators can be 
treated analytically. We estimate the common frequency 
for a particular case in Section III, and although similar 
analysis cannot be carried out for chaotic oscillators, this 
provides some insight into how the phenomenon arises 
more generally. The paper concludes with a discussion 
and summary in Section IV. 

II. AMPLITUDE DEATH 

The basic phenomenology of amplitude death in limit- 
cycle systems has been described in some detail Q,H|. We 
consider the general case of chaotic coupled systems, 

x = f x (x)+g x (y(t-r),x(t)) 

y = f y (y)+h y (x(t-r),y(t)) (1) 

where x and y denote the variables of the two subsys- 
tems. The dynamical equations are specified by f x and 
f y respectively, and g x and h y specify the couplings, r 
being the time-delay. 

Several scenarios are possible. In the simplest case, 
the subsystems specified by the variables x and y can be 
identical. This is dealt with in subsection A below. By 
altering the parameters, the two subsystems can be made 
nonidentical, and this case is treated in subsection B. Fi- 
nally, the two chaotic systems can be completely differ- 
ent, and this case is presented in subsection C. In all cases 
we demonstrate that amplitude death is possible, though 
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the mechanism in each case is somewhat different. For 
simplicity, the results are presented through applications 
to specific Rossler or Lorenz dynamical systems which 
have been extensively studied in the context of chaotic 
dynamics |12|. The present results are mainly numeri- 
cal since analytic results for such nonlinear systems are 
difficult to obtain. 

In each of the cases, we address the possibility of am- 
plitude suppression, namely the stabilization of limit cy- 
cles. For this phenomenon, some analysis is possible, and 
is presented in Sec. IIIII 



A. Identical Chaotic Subsystems 

The Rossler system ^| is a simple mathematical 
model of chemical kinetics that incorporates reaction- 
diffusion, and has been extensively studied in the past 
two decades as one of the simplest chaotic flows. Con- 
sider two identical Rossler oscillators denoted by vari- 
ables x and y 
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which are symmetrically coupled through x 2 and y 2 . We 
take numerical values of the parameters as a = b = 0.1 
and Ci = c 2 = c = 14 with the coupling parameters e and 
the delay r being treated as variable. 

In absence of coupling, e = 0, the subsystems, which 
are chaotic |13| , evolve independently. For finite coupling 
strength various states of motion are observed. These 
are indicated in the schematic Fig. ^a) within a rep- 
resentative range of coupling parameters e and r. The 
region marked C which is shaded black, corresponds to 
chaotic states while the white region shows the regions of 
regular behavior, namely periodic (P), fixed point (FP) 
and hypertorus (HT ) [lij dynamics. Numerical details 
are given in |l5l 1 1 6| . Finer analysis at a fixed value of 
coupling strength, e — 0.5 is shown in Fig. ^b). The 
three largest Lyapunov exponents are shown as function 
of the time delay r. The dotted line in Fig.^a) indicates 
the transition where the largest Lyapunov exponent (Ai) 
becomes negative [see Fig. ^b)]. The vertical arrow in 
Fig. ^b) denotes the parameter at which the slope of Ai 
changes sign. In Fig.^a), the locus of this point is plot- 
ted as a dashed line. The solid line in Fig. ^a) shows 



the transition from periodic solution to the hypertorus 
where Ai « A2 ~ 0. 

In the region marked P in Fig. ^b), the dynamics of 
the two subsystems settles onto limit cycles. Trajectories 
of the individual oscillators in the xi — x 2 plane with 
t = 0.6 are shown in Figs.^a) and (b) respectively. Here 
the trajectories of both the oscillators overlap since they 
are in complete synchronization. Further, these are also 
in phase since the phase difference is zero. There is also 
another periodic region between the FP and HT motions 
(between dotted and solid lines in Fig. ^a)) where tra- 
jectories of individual oscillators at r = 2.25 are shown in 
Figs.[2Ic) and (d). In the latter case both the oscillators 
show the same periodic motion (frequency locked states) 
(see Fig.|2Ic)) but they are out of phase [Fig.^d)], with 
a constant phase difference of ir. 

Both types of periodic motions (either in or out of 
phase) can arise from the following reasons 0|- Either 
the coupling stabilizes one of the (infinitely many) un- 
stable periodic orbits (UPO) 0] which are embedded in 
the chaotic set, or creates new periodic solutions with 
time period (T) different from the time delay, i.e. r =/= T 
(the perturbation term does not vanish here). In the for- 
mer case only those UPOs will be stabilized whose time 
periods are equal to the time delay, namely r = T . A 
trajectory will settle on an UPO and remain there since 
the feedback function becomes zero: the coupling term 
vanishes when y 2 (t — r) becomes equal to x 2 (t) . This im- 
plies that the perturbation does not change the original 
solution. In the latter situation, the effect of time delay 
is shown in Figs. |3[a-d). 

As t is increased further (between the dotted lines in 
Fig. ^a), where the largest Lyapunov exponent, Ai is 
negative as can be seen in Fig. ^b)) fixed point (FP) 
solutions are obtained. Examples of such solutions, be- 
fore and after the marked arrow in Fig. ^b) are shown in 
Figs.|2fe,f) and (g,h) respectively for t = 1.5 and r = 2. 
In both cases transient trajectories start spiraling into 
one of the fixed points xi+ = t/u = ax 3 +, x 2 * — y 2i , = 
-x 3 +,x 3 ± = y 3 *,x 3 + = [c- v/(c 2 - Aab)]/2a. In these 
regions, since the subsystems are identical, the coupling 
term vanishes and the fixed points remain the same as 
that of the uncoupled systems. The largest Lyapunov ex- 
ponent for the fixed point solution (Afp), the position of 
which does not change with time delay at a given value 
of the coupling strength, is also shown in Fig. ^b) (red 
line with circles). Since this is negative in the FP region 
(where Ai = A2 = Xfp), the fixed point is stable. Any 
attempt by the subsystems to move from this state in- 
creases the coupling term and brings the entire system 
back to the fixed point. Thus by using time delay cou- 
pling, amplitude death phenomena can be observed even 
in chaotic systems in a manner which is very similar to 
that in limit-cycles 0, • 

In the transition from chaos to amplitude death that 
occurs via a limit cycle, the mechanism is the same as 
that of the coupled limit cycle case (see also Sec. 
IIIip i.e. the transition from limit cycle to FP takes place 
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because a pair of complex eigenvalues cross the imaginary 
axis from right to left and the Lyapunov exponent Xfp 
becomes negative. Simultaneously the stability of the 
fixed point is lost as a pair of eigenvalues cross the axis 
from left to right. Thus, in this case amplitude death is 
initiated at a Hopf bifurcation. 

There are however some differences in the manner in 
which the fixed points are reached. The largest Lyapunov 
exponent has a change in slope at the point r c marked 
by the arrow, and differences can be seen in the phase 
relationship of the two oscillators: for r < r c [inset of 
Fig. Hff)] the two are in-phase, while for r > r c [inset 
of Fig. IHh)], they are out of phase. Shown in Fig. Qfc) 
is the phase difference (A9) between oscillators which is 
defined as A9 = (\9i(t) — 02(t)|) where (•} denotes the 
average over time while 9 At) ~ tan~ 1 [(x 2 (t) /xi(t)] and 
9 2 {t) ~ tan -1 [y 2 (t)/i/i(f)] HE El This clearly indicates 
that below r c the phase difference is zero while after it is 
7T. We also observe (see Figs.Etb,d,f,h)) that the two os- 
cillators are frequency-locked. Their common frequency, 
fl, (measured from the peak-to-peak separation |2J|) is 
shown in Fig. dd) where a step-like change is observed 
at r c . This behaviour is analyzed in a simpler system in 
the following Section. 

The change in the phase difference and common fre- 
quency occurs at parameter values corresponding to max- 
imal stability, namely when the Lyapunov exponent is at 
a local minimum. While details of the precise mechanism 
in chaotic systems needs further analysis, it can be easily 
verified that this behavior is quite general, and a similar 
transition occurs in a variety of model systems. 

Consider coupled identical chaotic Lorenz oscillators 
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with a = 10, ri = r 2 = 28, d = 8/3. The phase diagram 
as a function of the coupling parameters r and e is shown 
in Fig.|3Ja), while the Lyapunov exponents with time de- 
lay r at fixed e = 0.5 are shown in Fig. E^b). The step 
transition in phase difference and common frequency are 
shown in Fig. Etc) and (d) respectively. The behaviour is 
very similar to the Rossler case, Fig.^ Transient trajec- 
tories that go to the fixed point across the marked vertical 
arrow at r = 0.12 and r = 0.17 are shown in Figs.^e) 
(in-phase) and Figs.^f) (out of phase) respectively. The 
fixed point at (±y/d(r — 1), ±a/ 'd(r — 1), r — 1) is stable 



in the region marked FP, and has a negative Lyapunov 
exponent, Xfp 24]. In this region X±, = \ 2 = Xfp- Thus 
in this regime amplitude death also occurs. 

However there is difference in manner of transition 
from that of Fig. ^ I n the case of coupled limit-cycle 
P| or Rossler system, Eq. (J2J, transition is from peri- 
odic to fixed point while in this case it happens directly 
from chaotic dynamics to a fixed point. This implies 
that although neither UPOs are stabilized, nor new peri- 
odic solutions are created via the time-delay interaction, 
the fixed point is stabilized directly and this leads to an 
abrupt change in the largest Lyapunov exponent of the 
system, Ai. Such a transition, which has not been dis- 
cussed earlier, suggests that a new mechanism may be 
operative in amplitude death. 



B. Nonidentical Chaotic Subsystems 

We consider the effect of a difference in parameters be- 
tween the two subsystems. This is particularly important 
with respect to experimental realization of such phenom- 
ena since in practice it is impossible to ensure that all 
parameters of two different systems are exactly equal. 
Introduce a mismatch in the parameters of the Rossler 
system, Eq. (J2J. For c\ — 14 and c 2 — 18, various possi- 
ble states are shown in Fig. Ua). Other parameters are 
the same as in Fig. ^ A major difference, compared 
with Fig. d is the absence of hypertorus dynamics. The 
spectrum of Lyapunov exponents at e = 0.5 are shown 
in Fig. EJb); in the fixed point region all the Lyapunov 
exponents are negative. The periodic states at r = 0.25 
and r = 2.5 are shown in Figs. EJe) and (f) where the 
motions are in and out of phase respectively (these are in 
the frequency locked regime where time periods are same 
for both the oscillators). 

Shown in Figs. 01(g) and (h) are transient trajectories 
at t = 1 and r = 2 respectively. The insets are enlarged 
views of the corresponding transient and fixed point mo- 
tions. Here coupling introduces a new set of fixed points 
and the system settles in one of these. The fixed points 
for Eq. (J2J simply turn out to be 

xu = ax 3i , - e(x 3 * - y 3it ), 

x 2 ± = -x 3 +, 

x 3 ± = [(e-a)y3* + c 22/3*-&]/«/3*, 

2/1* = ay 3 + - e(y 3 ± -x 3it ), 

V2± = -2/3*, (4) 

where y 3i , is the root of the polynomial 

b + 2/3*[a2/3* - AV3+ - £3*) - c 2 ] = 0. 

The largest Lyapunov exponent Xfp (red line with 
circles), of the fixed point Xx* = 0.00135.., x 2 * = 
-0.0056.., x 3 ± = 0.0056.., y u = -7.914.. x lO" 5 ,^ = 
-0.00714.., and y 3 * = 0.00714.. is shown in Fig. jjb). 
Here the transition is from chaos to a fixed point occurs 
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indirectly, namely via a limit cycle at a Hopf bifurca- 
tion as in the case of identical coupled Rossler oscillators, 

Fig.n 

Similar step transitions in the phase difference and 
common frequency as in Figs. IHc,d) and|3](c,d) are also 
observed here in Figs. EJc,d). However the phase differ- 
ence are only « and ss it before and after the transition 
(marked by the arrow). 

Amplitude death phenomena thus also occurs in non- 
identical chaotic subsystems, though there can be subtle 
differences in the nature of the final state which can have 
either in-phase 0, or out of phase (2(| similar to the 
identical subsystem case. Similar results have been ob- 
tained [25] for the case of mismatched coupled Lorenz 
systems. 

C. Mixed Chaotic Subsystems 



Regions, denoted by R, where wild fluctuations (see 
Fig. Eta)) in the Lyapunov exponent occur are indica- 
tive of the presence of coexisting attractors with compli- 
cated basins [2{J |2]j ■ I have verified the presence of rid- 
dled [28| basins (results are not presented here) in this 
region, and identify at least two coexisting attractors: 
the chaotic trajectories (shown separately for Rossler and 
Lorenz subsystems in Figs. Etg) and (h) respectively), 
and the fixed point attractor, in Figs.Et e ) and (f). Here 
the transition is directly from chaos to fixed point similar 
to Lorenz oscillator, Eq. © (neither is an UPO stabliscd 
nor is a new periodic solution created) but in the present 
case the FP solution is created via riddling. 

III. AMPLITUDE SUPPRESSION: 
ANALYTICAL ESTIMATION OF THE 
FREQUENCY JUMP 



Since amplitude death occurs for interacting oscillators 
with very different time-periods, it is natural to consider 
the situation when the two subsystems are dynamically 
distinct. This is of even greater relevance experimentally, 
since a variety of different oscillator systems can be mu- 
tually coupled in a given natural situation. 

To examine the behavior that obtains when distinct os- 
cillators are coupled, I study coupled Rossler and Lorenz 
oscillators 
dxi{t) 
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(5) 



where parameters c\ = 18 and r 2 — 28. Other parame- 
ters are the same as in Eqs. J2J and A typical phase 
diagram is shown in Fig. Eta), and the Lyapunov expo- 
nent as a function of delay r is shown in Fig.Etb). Vari- 
ous possible states are indicated: the region of fixed point 
motion, i.e. amplitude death phenomena, is observed in 
the ranger ~ 0.3— r ~ 0.52 (where Ai = Xfp)- Trajecto- 
ries for r = 0.35 and r = 0.45 in amplitude death region 
are shown inEJe) and (f) respectively. The inset figures 
indicate that both the subsystems are neither completely 
in nor completely out of phase i.e. phase differences are 
away from zero or near ir (see Figs. Etc) and (d)) due to 
the completely distinct natures of the chaotic dynamics of 
the individual subsystems. Even though the subsystems 
are distinct, it has been possible to observe indications 
of a vestige of the transition in the phase difference when 
Lyapunov exponent is at a local minimum |25j. 



The phenomenon of frequency jump which can be ob- 
served (See Sections IIA, B) in identical Rossler and 
Lorenz chaotic systems can be seen in coupled limit cy- 
cles systems as well. Consider the case that has been 
studied in detail by Reddy et al. 0, 

Zx{t) = (l+iwi-|Zi(t)| 2 )Zi(t) + e[Z 2 (t-T)-Zi(t)] 



Z 2 (t) = {l + iu 2 -\Z 2 {i)\ 2 )Z 2 {t) + e[Z x {t-T) 



(6) 



where Zj {t),j = 1,2 are the complex amplitudes of j—th. 
oscillator with \Zj\ — 1 and corresponding angular fre- 
quencies io j. e is the coupling strength. For simplicity 
we consider identical oscillators with io\ — uj 2 . 

Specific results are shown for the case of lj± = 9, which 
was the case used in Ref. p. The origin is a fixed point, 
i.e. Zj = 0. The spectrum of Lyapunov exponents and 
A fp the exponent (the red line with circles) for the fixed 
point, for coupling strength e = 10 is shown in Fig. Eta). 
The resulting limit cycles at r = 0.05 are in phase, and 
at t = 0.25 are out of phase, and these are shown in 
Fig. Etb) and (c) respectively. The amplitude death re- 
gion is in between r ~ 0.1 and 0.2 where all Lyapunov ex- 
ponents are negative, Xi = X 2 = Xfp- The transient be- 
haviour of in and out of phase motions in this amplitude 
death region are shown in Fig. El (d) and (e) at r = 0.13 
and 0.19 respectively. The transition form limit cycle to 
fixed point motion (say from r = 0.05 to r = 0.14) is 
happening at the Hopf bifurcation. The Lyapunov ex- 
ponent of the fixed point, Xfp, indicates this transition 
clearly, when it becomes negative. 

Similar to Figs. Qtc) and Et c )) shown in Fig. Etf) 
is the phase difference between oscillators of Eq. ©■ 
The phases of the individual oscillators are defined here 
as 6j = taa [Im(Zj) / Re(Zj)]. The phase difference 
changes drastically from to ir when the largest Lya- 
punov exponent is minimum (the system is in most 
stable stat). 
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There is also a similar step transition in frequency 
across which in- and out-of — phase motions are present. 
The numerically calculated frequency, Q, of Eq. JBJ with 
time delay, is shown by the dashed line in Fig. Etg) • In 
order to understand this transition I consider the charac- 
teristic eigenvalue equation for Eq. JfJJl, by assuming the 
linear perturbation to vary as e A *, which can be written 
as 0, 

A 2 - 2(a + iuj)X + (a 2 - uj 2 + i2auj) - e 2 e" 2Ar = 0, 

(7) 

where a = 1 — e. 

Letting A = a + i[3 where a and (3 are real and imag- 
inary part of the eigenvalues in Eq. Q, and separating 
real and imaginary parts leads to the equations 

a 2 - f3 2 - 2{aa - 0u) + a 2 - uj 2 - e 2 e~ 2aT cos (2/3t) = 

(8) 

and 

2a(3-2(au + afi) + 2auj + e 2 e- 2aT s\n{2(5T) = 0. 

(9) 

As these equations are difficult to solve analytically 
we use the Lyapunov exponent of the fixed point, Xfp, 
to find (3 numerically. The Lyapunov exponent of a 
fixed point is equal to the real part of the eigenvalue, 
i.e. Xfp — a, we use this in Eqs. 1® and ©■ The re- 
sulting roots of these equations give (5 and are shown 
separately by circles and stars in Fig. Elg) ■ The final so- 
lution which satisfying both these equations will be the 
solution of the eigenvalue equation, Eq. Q) . Here we see 
that prior to the transition, the smaller root of Eq. © 
matches Eq. © while after the transition, it is the larger 
root. The agreement with (3 computed directly from the 
eigenvalue equation and numerically calculated frequency 
(dashed line) is excellent. 

This analysis has been possible for the case of coupled 
limit cycle oscillators. Although it cannot be extended 
to the case of chaotic systems, it is possible to conjecture 
that the underlying mechanism could be similar, namely 
that the transition is due to a jump in the imaginary part 
of the eigenvalue of the fixed point. 

IV. DISCUSSION AND SUMMARY 

In this paper I have shown that the phenomenology 
developed for the cause of amplitude death in systems 
with limit-cycle oscillators jlj, |2j, |3| can be extended to 
the case of chaotic dynamical systems. This also ex- 
tends earlier work on amplitude modulation [29|, where 
time-delay coupling has been used to stabilize undesir- 
able low-frequency chaotic fluctuations in a semiconduc- 
tor laser. (There it was seen there that time-delay cou- 
pling could convert chaotic oscillations of the laser field 



with power dropouts into quasiperiodic motion without 
power dropouts.) 

Amplitude death in chaotic systems is a consequence 
of time-delays that stabilize fixed points that might ei- 
ther have existed in the uncoupled system, or that might 
have been created through the interaction. The evidence 
presented here is largely numerical. We considered sepa- 
rately the cases of the interacting systems being (a) iden- 
tical, (b) the same but with mismatched parameters, and 
(c) completely distinct. In all cases, amplitude death 
can occur, and moreover, this can occur over a range of 
parameter values. The mechanics through which chaos 
converts into fixed point motion are different in these 
example: some cases it happens via limit-cycle at Hopf 
bifurcation while other cases fixed point gets stablised 
directly. In mixed systems this stabilization is initiated 
via riddling phenomenon. 

We have used the Lyapunov exponent as indicator to 
detect the different dynamical regimes, and to detect 
transitions, as for example, in the case of identical os- 
cillators, when the subsystems are oscillating in or out of 
phase with each other. When the subsystems are max- 
imally stable, namely the largest Lyapunov exponent is 
minimum, a "phase" transition occurs, and motion goes 
from being in phase to being out of phase. Measures 
that detect this transition are the phase difference and 
frequency; although not possible to analyze in chaotic 
systems, it is possible to get some understanding of this 
transition for a limit-cycle case. Apart from its mathe- 
matical interest, the possibility of using a transition from 
dynamics on a limit cycle which is in phase to that which 
is out of phase in practical applications suggest itself. In 
experimental application of amplitude death phenomena, 
where transient dynamics always persists, it will be im- 
portant to know the parameters where motion is cither 
in-phase (low frequency) or out of phase (high frequency). 
An im por tant possible use of it could be a coupled chaotic 
lasers [2!j , where relatively higher degree of constant out- 
put could be obtained by keeping the parameters of the 
lasers and the time-delay coupling in the out of phase 
regime compared to that of the in-phase [2{| . 

Experimental verification of this generalization of the 
conventional amplitude-death phenomenon that occurs 
in limit-cycle oscillators may in fact, prove to be techni- 
cally simple to achieve since it persists even after param- 
eters mismatch, and occurs for different types of coupled 
oscillators, bypassing the need for irregular oscillations 
(e. g. Ref. [2|| ). 
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FIG. 1: (Color online) (a) Schematic phase diagram of iden- 
tical Rossler oscillators, Eq. (|5J, in the e — r plane for 
ci = C2 = 18. In the dark region the motion is chaotic (C) 
while in the white region there are regular states of different 
kind e. g. fixed point (FP), periodic cycles (P) and hyper- 
torus (HT). Numerical details are given in |15|. The dotted 
line indicates the locus where largest Lyapunov exponent (Ai) 
becomes negative, (b) The spectrum of Lyapunov exponents 
(black, green and blue correspond to Ai, A2, and A3 respec- 
tively) as a function of time delay r at fixed coupling strength 
e — 0.5. The vertical arrow in (a) shows the parameter value, 
T c , when the motion goes from being in phase to out of phase. 
The largest Lyapunov exponent of the fixed point is shown in 
red. In the amplitude death range all the Lyapunov expo- 
nents are negative and Ai = A2 = Xfp- The transition is 
from limit cycle to a fixed point, (c) The phase difference 
between oscillators with time delay r, which before and after 
the transition are equal to and n. (d) Numerically calcu- 
lated common frequency, f2, as a function of the time delay 
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FIG. 2: Trajectories of oscillator 1 (solid line) and 2 (dashed 
line) in phase space (left panel) and with time (right panel) 
(a) and (b) at t — 0.6, (c) and (d) at r = 2.25, (d) and (e) at 
r = 1.5, and, (g) and (h) at r = 2 respectively. 
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FIG. 3: (Color online) (a) Schematic phase diagram for the 
Lorenz system Eq. JJ) in e — r plane, (b) Spectrum of Lya- 
punov exponents (black, green and blue corresponding to Ai , 
A2, and A3 respectively) as a function of the time delay, r, 
at fixed coupling strength e = 0.5. The other details remains 
the same as that in the Fig. In the amplitude death range 
all the Lyapunov exponents are negative and Ai = A2 = Xfp', 
the transition is directly from chaos to a fixed point, (c) The 
phase difference between oscillators with time delay r. The 
phase difference before and after the transition are equal to 
and 7T. (d) Numerically calculated common frequency, fi. 
Trajectories of oscillators 1 fsolid line"! and 2 ("dashed line"! 
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FIG. 4: (Color online) (a) Schematic phase diagram of Eq. J2J 
in e — t plane. Parameters values are ci = 14 and C2 = 18. 
(b) Spectrum of Lyapunov exponent (black, green and blue 
correspond to Ai, A2, and A3 respectively) as a function of the 
time delay, r, at fixed coupling strength e — 0.5. Other details 
remains the same as that of Fig. Q In amplitude death range 
all Lyapunov exponents are negative and Ai = A2 = \fp- 
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FIG. 5: (Color online) (a) Schematic phase diagram for the 
mixed Rossler and Lorenz system, Eq. (J^J, in e — r plane, 
(b) The largest Lyapunov exponent (in black) as a function 
of time delay, r, at fixed coupling strength e = 0.5. R de- 
notes regions of wild fluctuations in the Lyapunov exponent 
where chaotic and stable fixed point solutions coexist. The 
red line with circles represents the largest Lyapunov exponent 
of a fixed point (xx* = 4.200308.., x 2 * = -0.007246.., x 3ir = 
-x 2 +,yi+ = -8.406415.., y 2ir = yi*, and y 3 * = 26.500431..). 
In contrast to previous examples, there is no discontinuous 
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FIG. 6: (Color online) (a) Spectrum of Lyapunov exponents 
for the coupled limit cycle oscillators (black and green corre- 
spond to Ai and A2 respectively) with time delay, r, at fixed 
coupling strength e — 10. The red line with circles is the 
largest Lyapunov exponent for the fixed point, Zj — 0. In 
amplitude death region all Lyapunov exponents are negative 
and Ai = A2 = \fp- Trajectories of the oscillators 1 (solid 
line) and 2 (dashed line) with limit cycles at (b) r = 0.05 
(in phase) and (c) r = .25 (out of phase). ReZj represents 
the real part for Zj of j—th oscillator. Transient trajectories 
of oscillators 1 (solid line) and 2 (dashed line) with time at 
(d) r = 0.14 (in phase) and (e) r = 0.19 (out of Phase) re- 
spectively for fixed point solutions, (f) The phase difference 
between oscillators with time delay r which is equal to and 
7r before and after the transition respectively, (g) Solution of 
Eqs. © (black circles) and © (red stars) with parameter r. 
Blue dashed line indicates the numerically calculated common 
frequency, fl. 



